Gramian matrix

In linear algebra, the Gramian matrix (or Gram matrix or Gramian) of a set of vectors $v_1,\dots, v_n$ in an inner product space is the Hermitian matrix of inner products, whose entries are given by $G_{ij}=\langle v_j, v_i \rangle$ .

An important application is to compute linear independence: a set of vectors is linearly independent if and only if the Gram determinant (the determinant of the Gram matrix) is non-zero.

It is named after Jørgen Pedersen Gram.

1 Examples
- 1.1 Applications
2 Properties
- 2.1 Positive semidefinite
- 2.2 Change of basis
3 Gram determinant
4 References
5 External links

Examples

Most commonly, the vectors are elements of an Euclidean space, or are functions in an L² space, such as continuous functions on a compact interval [a, b] (which are a subspace of L²([a, b])).

Given real-valued functions $\{\ell_i(\cdot),\,i=1,\dots,n\}$ on the interval $[t_0,t_f]$ , the Gram matrix $G=[G_{ij}]$ , is given by the standard inner product on functions:

$G_{ij}=\int_{t_0}^{t_f} \ell_i(\tau)\ell_j(\tau)\, d\tau.$

Given a real matrix A, the matrix A^TA is a Gram matrix (of the columns of A), while the matrix AA^T is the Gram matrix of the rows of A.

For a general bilinear form B on a finite-dimensional vector space over any field we can define a Gram matrix G attached to a set of vectors $v_1,\dots, v_n$ by $G_{i,j} = B(v_i,v_j) \,$ . The matrix will be symmetric if the bilinear form B is symmetric.

Applications

If the vectors are centered random variables, the Gramian is proportional to the covariance matrix, with the scaling determined by the number of elements in the vector.
In quantum chemistry, the Gram matrix of a set of basis vectors is the overlap matrix.
In control theory (or more generally systems theory), the controllability Gramian and observability Gramian determine properties of a linear system.
Gramian matrices arise in covariance structure model fitting (see e.g., Jamshidian and Bentler, 1993, Applied Psychological Measurement, Volume 18, pp. 79–94).
In the finite element method, the Gram matrix arises from approximating a function from a finite dimensional space; the Gram matrix entries are then the inner products of the basis functions of the finite dimensional subspace.
In machine learning, kernel functions are often represented as Gram matrices.^[1]

Properties

Positive semidefinite

The Gramian matrix is positive semidefinite, and every positive semidefinite matrix is the Gramian matrix for some set of vectors. Further, in finite-dimensions it determines the vectors up to isomorphism, i.e. any two sets of vectors with the same Gramian matrix must be related by a single unitary matrix. The Gramian matrix of any orthonormal basis is the identity matrix.

The infinite-dimensional analog of this statement is Mercer's theorem.

Change of basis

Under change of basis represented by an invertible matrix P, the Gram matrix will change by a matrix congruence to P^TGP.

Gram determinant

The Gram determinant or Gramian is the determinant of the Gram matrix:

$G(x_1,\dots, x_n)=\begin{vmatrix} \langle x_1,x_1\rangle & \langle x_1,x_2\rangle &\dots & \langle x_1,x_n\rangle\\ \langle x_2,x_1\rangle & \langle x_2,x_2\rangle &\dots & \langle x_2,x_n\rangle\\ \vdots&\vdots&\ddots&\vdots\\ \langle x_n,x_1\rangle & \langle x_n,x_2\rangle &\dots & \langle x_n,x_n\rangle\end{vmatrix}.$

Geometrically, the Gram determinant is the square of the volume of the parallelotope formed by the vectors. In particular, the vectors are linearly independent if and only if the Gram determinant is nonzero (if and only if the Gram matrix is nonsingular).

The Gram determinant can also be expressed in terms of the exterior product of vectors by

$G(x_1,\dots,x_n) = \| x_1\wedge\cdots\wedge x_n\|^2.$

References

^ Lanckriet, G. R. G., N. Cristianini, P. Bartlett, L. E. Ghaoui, and M. I. Jordan. “Learning the kernel matrix with semidefinite programming.” The Journal of Machine Learning Research 5 (2004): 29.

Barth, Nils (1999). "The Gramian and K-Volume in N-Space: Some Classical Results in Linear Algebra". Journal of Young Investigators 2. http://www.jyi.org/volumes/volume2/issue1/articles/barth.html.

External links

Volumes of parallelograms by Frank Jones